function U = upsilon(q, param)

    sigma = param.sigma;
    epsilon = param.epsilon;
        
   U = 1 + (sigma-1)*exp(1/epsilon)*epsilon^(sigma/epsilon-1)*...
        (uig(sigma/epsilon,1/epsilon) - ...
        uig(sigma/epsilon, (q.^(epsilon/sigma))/epsilon));
    
end
